Late Pleistocene stickleback environmental genomes reveal the chronology of freshwater adaptation.

Directly observing the chronology and tempo of adaptation in response to ecological change is rarely possible in natural ecosystems. Sedimentary ancient DNA (sedaDNA) has been shown to be a tractable source of genome-scale data of long-dead organisms1,2,3 and to thereby potentially provide an understanding of the evolutionary histories of past populations.4,5 To date, time series of ecosystem biodiversity have been reconstructed from sedaDNA, typically using DNA metabarcoding or shotgun sequence data generated from less than 1 g of sediment.6,7 Here, we maximize sequence coverage by extracting DNA from ∼50× more sediment per sample than the majority of previous studies1,2,3 to achieve genotype resolution. From a time series of Late Pleistocene sediments spanning from a marine to freshwater ecosystem, we compare adaptive genotypes reconstructed from the environmental genomes of three-spined stickleback at key time points of this transition. We find a staggered temporal dynamic in which freshwater alleles at known loci of large effect in marine-freshwater divergence of three-spined stickleback (e.g., EDA)8 were already established during the brackish phase of the formation of the isolation basin. However, marine alleles were still detected across the majority of marine-freshwater divergence-associated loci, even after the complete isolation of the lake from marine ingression. Our retrospective approach to studying adaptation from environmental genomes of three-spined sticklebacks at the end of the last glacial period complements contemporary experimental approaches9,10,11 and highlights the untapped potential for retrospective "evolve and resequence" natural experiments using sedaDNA.

Reproductive isolation in a three-way contact zone.

Contact zones between divergent forms within a species provide insight into the role of gene flow in adaptation and speciation. Previous work has focused on contact zones involving only two divergent forms, but in nature, many more than two populations may overlap simultaneously and experience gene flow. Patterns of introgression in wild populations are, therefore, likely much more complicated than is often assumed. We begin to address this gap in current knowledge by investigating patterns of divergence and introgression across a complex natural contact zone. We use phenotypic and genomic data to confirm the existence of a three-way contact zone among divergent freshwater resident, saltwater resident and saltwater migratory three-spined stickleback (Gasterosteus aculeatus) on the island of North Uist, Scottish Western Isles. We find evidence for hybridization, mostly between saltwater resident and saltwater migratory forms. Despite hybridization, genomic analyses reveal pairwise islands of divergence between all forms that are maintained across the contact zone. Genomic cline analyses also provide evidence for selection and/or hybrid incompatibilities in divergent regions. Divergent genomic regions occur across multiple chromosomes and involve many known adaptive loci and several chromosomal inversions. We also identify distinct immune gene expression profiles between forms, but no evidence for transgressive expression in hybrids. Our results suggest that reproductive isolation is maintained in this three-way contact zone, despite some hybridization, and that reduced recombination in chromosomal inversions may play an important role in maintaining this isolation.

Genomic changes underlying repeated niche shifts in an adaptive radiation.

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space.

Longer or shorter spines: Reciprocal trait evolution in stickleback via triallelic regulatory changes in Stanniocalcin2a.

Vertebrates have repeatedly modified skeletal structures to adapt to their environments. The threespine stickleback is an excellent system for studying skeletal modifications, as different wild populations have either increased or decreased the lengths of their prominent dorsal and pelvic spines in different freshwater environments. Here we identify a regulatory locus that has a major morphological effect on the length of stickleback dorsal and pelvic spines, which we term Maser (major spine enhancer). Maser maps in a closely linked supergene complex that controls multiple armor, feeding, and behavioral traits on chromosome IV. Natural alleles in Maser are differentiated between marine and freshwater sticklebacks; however, alleles found among freshwater populations are also differentiated, with distinct alleles found in short- and long-spined freshwater populations. The distinct freshwater alleles either increase or decrease expression of the bone growth inhibitor gene Stanniocalcin2a in developing spines, providing a simple genetic mechanism for either increasing or decreasing spine lengths in natural populations. Genomic surveys suggest many recurrently differentiated loci in sticklebacks are similarly specialized into three or more distinct alleles, providing multiple ancient standing variants in particular genes that may contribute to a range of phenotypes in different environments.

Predicting future from past: The genomic basis of recurrent and rapid stickleback evolution.

Similar forms often evolve repeatedly in nature, raising long-standing questions about the underlying mechanisms. Here, we use repeated evolution in stickleback to identify a large set of genomic loci that change recurrently during colonization of freshwater habitats by marine fish. The same loci used repeatedly in extant populations also show rapid allele frequency changes when new freshwater populations are experimentally established from marine ancestors. Marked genotypic and phenotypic changes arise within 5 years, facilitated by standing genetic variation and linkage between adaptive regions. Both the speed and location of changes can be predicted using empirical observations of recurrence in natural populations or fundamental genomic features like allelic age, recombination rates, density of divergent loci, and overlap with mapped traits. A composite model trained on these stickleback features can also predict the location of key evolutionary loci in Darwin's finches, suggesting that similar features are important for evolution across diverse taxa.

Ancient and modern stickleback genomes reveal the demographic constraints on adaptation.

Adaptation is typically studied by comparing modern populations with contrasting environments. Individuals persisting in the ancestral habitat are typically used to represent the ancestral founding population; however, it has been questioned whether these individuals are good proxies for the actual ancestors.1 To address this, we applied a paleogenomics approach2 to directly access the ancestral genepool: partially sequencing the genomes of two 11- to 13,000-year-old stickleback recovered from the transitionary layer between marine and freshwater sediments of two Norwegian isolation lakes3 and comparing them with 30 modern stickleback genomes from the same lakes and adjacent marine fjord, in addition to a global dataset of 20 genomes.4 The ancient stickleback shared genome-wide ancestry with the modern fjord population, whereas modern lake populations have lost substantial ancestral variation following founder effects, and subsequent drift and selection. Freshwater-adaptive alleles found in one ancient stickleback genome have not risen to high frequency in the present-day population from the same lake. Comparison to the global dataset suggested incomplete adaptation to freshwater in our modern lake populations. Our findings reveal the impact of population bottlenecks in constraining adaptation due to reduced efficacy of selection on standing variation present in founder populations.

Fitness maps to a large-effect locus in introduced stickleback populations.

Mutations of small effect underlie most adaptation to new environments, but beneficial variants with large fitness effects are expected to contribute under certain conditions. Genes and genomic regions having large effects on phenotypic differences between populations are known from numerous taxa, but fitness effect sizes have rarely been estimated. We mapped fitness over a generation in an F2 intercross between a marine and a lake stickleback population introduced to a freshwater pond. A quantitative trait locus map of the number of surviving offspring per F2 female detected a single, large-effect locus near Ectodysplasin (Eda), a gene having an ancient freshwater allele causing reduced bony armor and other changes. F2 females homozygous for the freshwater allele had twice the number of surviving offspring as homozygotes for the marine allele, producing a large selection coefficient, s = 0.50 ± 0.09 SE. Correspondingly, the frequency of the freshwater allele increased from 0.50 in F2 mothers to 0.58 in surviving offspring. We compare these results to allele frequency changes at the Eda gene in an Alaskan lake population colonized by marine stickleback in the 1980s. The frequency of the freshwater Eda allele rose steadily over multiple generations and reached 95% within 20 y, yielding a similar estimate of selection, s = 0.49 ± 0.05, but a different degree of dominance. These findings are consistent with other studies suggesting strong selection on this gene (and/or linked genes) in fresh water. Selection on ancient genetic variants carried by colonizing ancestors is likely to increase the prevalence of large-effect fitness variants in adaptive evolution.

Genome-wide recombination map construction from single individuals using linked-read sequencing.

Meiotic recombination rates vary across the genome, often involving localized crossover hotspots and coldspots. Studying the molecular basis and mechanisms underlying this variation has been challenging due to the high cost and effort required to construct individualized genome-wide maps of recombination crossovers. Here we introduce a new method, called ReMIX, to detect crossovers from gamete DNA of a single individual using Illumina sequencing of 10X Genomics linked-read libraries. ReMIX reconstructs haplotypes and identifies the valuable rare molecules spanning crossover breakpoints, allowing quantification of the genomic location and intensity of meiotic recombination. Using a single mouse and stickleback fish, we demonstrate how ReMIX faithfully recovers recombination hotspots and landscapes that have previously been built using hundreds of offspring. ReMIX provides a high-resolution, high-throughput, and low-cost approach to quantify recombination variation across the genome, providing an exciting opportunity to study recombination among multiple individuals in diverse organisms.

Predominance of cis-regulatory changes in parallel expression divergence of sticklebacks.

Regulation of gene expression is thought to play a major role in adaptation, but the relative importance of cis- and trans- regulatory mechanisms in the early stages of adaptive divergence is unclear. Using RNAseq of threespine stickleback fish gill tissue from four independent marine-freshwater ecotype pairs and their F1 hybrids, we show that cis-acting (allele-specific) regulation consistently predominates gene expression divergence. Genes showing parallel marine-freshwater expression divergence are found near to adaptive genomic regions, show signatures of natural selection around their transcription start sites and are enriched for cis-regulatory control. For genes with parallel increased expression among freshwater fish, the quantitative degree of cis- and trans-regulation is also highly correlated across populations, suggesting a shared genetic basis. Compared to other forms of regulation, cis-regulation tends to show greater additivity and stability across different genetic and environmental contexts, making it a fertile substrate for the early stages of adaptive evolution.

Experimental evidence for rapid genomic adaptation to a new niche in an adaptive radiation.

A substantial part of biodiversity is thought to have arisen from adaptive radiations in which one lineage rapidly diversified into multiple lineages specialized to many different niches. However, selection and drift reduce genetic variation during adaptation to new niches and may thus prevent or slow down further niche shifts. We tested whether rapid adaptation is still possible from a highly derived ecotype in the adaptive radiation of threespine stickleback on the Haida Gwaii archipelago, Western Canada. In a 19-year selection experiment, we let giant sticklebacks from a large blackwater lake evolve in a small clearwater pond without vertebrate predators. A total of 56 whole genomes from the experiment and 26 natural populations revealed that adaptive genomic change was rapid in many small genomic regions and encompassed 75% of the change between 12,000-year-old ecotypes. Genomic change was as fast as phenotypic change in defence and trophic morphology, and both were largely parallel between the short-term selection experiment and long-term natural adaptive radiation. Our results show that functionally relevant standing genetic variation can persist in derived radiation members, allowing adaptive radiations to unfold very rapidly.

Detecting differential copy number variation between groups of samples.

We present a method to detect copy number variants (CNVs) that are differentially present between two groups of sequenced samples. We use a finite-state transducer where the emitted read depth is conditioned on the mappability and GC-content of all reads that occur at a given base position. In this model, the read depth within a region is a mixture of binomials, which in simulations matches the read depth more closely than the often-used negative binomial distribution. The method analyzes all samples simultaneously, preserving uncertainty as to the breakpoints and magnitude of CNVs present in an individual when it identifies CNVs differentially present between the two groups. We apply this method to identify CNVs that are recurrently associated with postglacial adaptation of marine threespine stickleback (Gasterosteus aculeatus) to freshwater. We identify 6664 regions of the stickleback genome, totaling 1.7 Mbp, which show consistent copy number differences between marine and freshwater populations. These deletions and duplications affect both protein-coding genes and cis-regulatory elements, including a noncoding intronic telencephalon enhancer of DCHS1 The functions of the genes near or included within the 6664 CNVs are enriched for immunity and muscle development, as well as head and limb morphology. Although freshwater stickleback have repeatedly evolved from marine populations, we show that freshwater stickleback also act as reservoirs for ancient ancestral sequences that are highly conserved among distantly related teleosts, but largely missing from marine stickleback due to recent selective sweeps in marine populations.

Genetic Coupling of Female Mate Choice with Polygenic Ecological Divergence Facilitates Stickleback Speciation.

Ecological speciation with gene flow is widespread in nature [1], but it presents a conundrum: how are associations between traits under divergent natural selection and traits that contribute to assortative mating maintained? Theoretical models suggest that genetic mechanisms inhibiting free recombination between loci underlying these two types of traits (hereafter, "genetic coupling") can facilitate speciation [2-4]. Here, we perform a direct test for genetic coupling by mapping both divergent traits and female mate choice in a classic model of ecological speciation: sympatric benthic and limnetic threespine stickleback (Gasterosteus aculeatus). By measuring mate choice in F2 hybrid females, we allowed for recombination between loci underlying assortative mating and those under divergent ecological selection. In semi-natural mating arenas in which females had access to both benthic and limnetic males, we found that F2 females mated with males similar to themselves in body size and shape. In addition, we found two quantitative trait loci (QTLs) associated with female mate choice that also predicted female morphology along the benthic-limnetic trait axis. Furthermore, a polygenic genetic model that explains adaptation to contrasting benthic and limnetic feeding niches [5] also predicted F2 female mate choice. Together, these results provide empirical evidence that genetic coupling of assortative mating with traits under divergent ecological selection helps maintain species in the face of gene flow, despite a polygenic basis for adaptation to divergent environments.

Convergent evolution of SWS2 opsin facilitates adaptive radiation of threespine stickleback into different light environments.

Repeated adaptation to a new environment often leads to convergent phenotypic changes whose underlying genetic mechanisms are rarely known. Here, we study adaptation of color vision in threespine stickleback during the repeated postglacial colonization of clearwater and blackwater lakes in the Haida Gwaii archipelago. We use whole genomes from 16 clearwater and 12 blackwater populations, and a selection experiment, in which stickleback were transplanted from a blackwater lake into an uninhabited clearwater pond and resampled after 19 y to test for selection on cone opsin genes. Patterns of haplotype homozygosity, genetic diversity, site frequency spectra, and allele-frequency change support a selective sweep centered on the adjacent blue- and red-light sensitive opsins SWS2 and LWS. The haplotype under selection carries seven amino acid changes in SWS2, including two changes known to cause a red-shift in light absorption, and is favored in blackwater lakes but disfavored in the clearwater habitat of the transplant population. Remarkably, the same red-shifting amino acid changes occurred after the duplication of SWS2 198 million years ago, in the ancestor of most spiny-rayed fish. Two distantly related fish species, bluefin killifish and black bream, express these old paralogs divergently in black- and clearwater habitats, while sticklebacks lost one paralog. Our study thus shows that convergent adaptation to the same environment can involve the same genetic changes on very different evolutionary time scales by reevolving lost mutations and reusing them repeatedly from standing genetic variation.

Evolving New Skeletal Traits by cis-Regulatory Changes in Bone Morphogenetic Proteins.

Changes in bone size and shape are defining features of many vertebrates. Here we use genetic crosses and comparative genomics to identify specific regulatory DNA alterations controlling skeletal evolution. Armor bone-size differences in sticklebacks map to a major effect locus overlapping BMP family member GDF6. Freshwater fish express more GDF6 due in part to a transposon insertion, and transgenic overexpression of GDF6 phenocopies evolutionary changes in armor-plate size. The human GDF6 locus also has undergone distinctive regulatory evolution, including complete loss of an enhancer that is otherwise highly conserved between chimps and other mammals. Functional tests show that the ancestral enhancer drives expression in hindlimbs but not forelimbs, in locations that have been specifically modified during the human transition to bipedalism. Both gain and loss of regulatory elements can localize BMP changes to specific anatomical locations, providing a flexible regulatory basis for evolving species-specific changes in skeletal form.

Extent of QTL Reuse During Repeated Phenotypic Divergence of Sympatric Threespine Stickleback.

How predictable is the genetic basis of phenotypic adaptation? Answering this question begins by estimating the repeatability of adaptation at the genetic level. Here, we provide a comprehensive estimate of the repeatability of the genetic basis of adaptive phenotypic evolution in a natural system. We used quantitative trait locus (QTL) mapping to discover genomic regions controlling a large number of morphological traits that have diverged in parallel between pairs of threespine stickleback (Gasterosteus aculeatus species complex) in Paxton and Priest lakes, British Columbia. We found that nearly half of QTL affected the same traits in the same direction in both species pairs. Another 40% influenced a parallel phenotypic trait in one lake but not the other. The remaining 10% of QTL had phenotypic effects in opposite directions in the two species pairs. Similarity in the proportional contributions of all QTL to parallel trait differences was about 0.4. Surprisingly, QTL reuse was unrelated to phenotypic effect size. Our results indicate that repeated use of the same genomic regions is a pervasive feature of parallel phenotypic adaptation, at least in sticklebacks. Identifying the causes of this pattern would aid prediction of the genetic basis of phenotypic evolution.

A recurrent regulatory change underlying altered expression and Wnt response of the stickleback armor plates gene EDA.

Armor plate changes in sticklebacks are a classic example of repeated adaptive evolution. Previous studies identified ectodysplasin (EDA) gene as the major locus controlling recurrent plate loss in freshwater fish, though the causative DNA alterations were not known. Here we show that freshwater EDA alleles have cis-acting regulatory changes that reduce expression in developing plates and spines. An identical T → G base pair change is found in EDA enhancers of divergent low-plated fish. Recreation of the T → G change in a marine enhancer strongly reduces expression in posterior armor plates. Bead implantation and cell culture experiments show that Wnt signaling strongly activates the marine EDA enhancer, and the freshwater T → G change reduces Wnt responsiveness. Thus parallel evolution of low-plated sticklebacks has occurred through a shared DNA regulatory change, which reduces the sensitivity of an EDA enhancer to Wnt signaling, and alters expression in developing armor plates while preserving expression in other tissues.

Genetics of ecological divergence during speciation.

Ecological differences often evolve early in speciation as divergent natural selection drives adaptation to distinct ecological niches, leading ultimately to reproductive isolation. Although this process is a major generator of biodiversity, its genetic basis is still poorly understood. Here we investigate the genetic architecture of niche differentiation in a sympatric species pair of threespine stickleback fish by mapping the environment-dependent effects of phenotypic traits on hybrid feeding and performance under semi-natural conditions. We show that multiple, unlinked loci act largely additively to determine position along the major niche axis separating these recently diverged species. We also find that functional mismatch between phenotypic traits reduces the growth of some stickleback hybrids beyond that expected from an intermediate phenotype, suggesting a role for epistasis between the underlying genes. This functional mismatch might lead to hybrid incompatibilities that are analogous to those underlying intrinsic reproductive isolation but depend on the ecological context.

Genomics and the origin of species.

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Phylogeography and adaptation genetics of stickleback from the Haida Gwaii archipelago revealed using genome-wide single nucleotide polymorphism genotyping.

Threespine stickleback populations are model systems for studying adaptive evolution and the underlying genetics. In lakes on the Haida Gwaii archipelago (off western Canada), stickleback have undergone a remarkable local radiation and show phenotypic diversity matching that seen throughout the species distribution. To provide a historical context for this radiation, we surveyed genetic variation at >1000 single nucleotide polymorphism (SNP) loci in stickleback from over 100 populations. SNPs included markers evenly distributed throughout genome and candidate SNPs tagging adaptive genomic regions. Based on evenly distributed SNPs, the phylogeographic pattern differs substantially from the disjunct pattern previously observed between two highly divergent mtDNA lineages. The SNP tree instead shows extensive within watershed population clustering and different watersheds separated by short branches deep in the tree. These data are consistent with separate colonizations of most watersheds, despite underlying genetic connections between some independent drainages. This supports previous suppositions that morphological diversity observed between watersheds has been shaped independently, with populations exhibiting complete loss of lateral plates and giant size each occurring in several distinct clades. Throughout the archipelago, we see repeated selection of SNPs tagging candidate freshwater adaptive variants at several genomic regions differentiated between marine-freshwater populations on a global scale (e.g. EDA, Na/K ATPase). In estuarine sites, both marine and freshwater allelic variants were commonly detected. We also found typically marine alleles present in a few freshwater lakes, especially those with completely plated morphology. These results provide a general model for postglacial colonization of freshwater habitat by sticklebacks and illustrate the tremendous potential of genome-wide SNP data sets hold for resolving patterns and processes underlying recent adaptive divergences.